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Abstract 

We present a new efficient analytical approximation scheme to two-point boundary 
value problems of ordinary differential equations (ODEs) adapted to the study of 
the derivative expansion of the exact renormalization group equations. It is based 
on a compactification of the complex plane of the independent variable using a 
mapping of an angular sector onto a unit disc. We explicitly treat, for the scalar 
field, the local potential approximations of the Wegner-Houghton equation in the 
dimension (i = 3 and of the Wilson- Polchinski equation for some values of d £ ]2, 3]. 
We then consider, for d = 3, the coupled ODEs obtained by Morris at the second 
order of the derivative expansion. In both cases the fixed points and the eigenvalues 
attached to them are estimated. Comparisons of the results obtained are made with 
the shooting method and with the other analytical methods available. The best 
accuracy is reached with our new method which presents also the advantage of 
being very fast. Thus, it is well adapted to the study of more complicated systems 
of equations. 



Key words: Exact renormalisation group. Derivative expansion, Critical 

exponents, Two-point boundary value problem 

PACS: 02.30.Hq, 02.30.Mv, 02.60.Lj, 05.10.Cc, ll.lO.Gh, 64.60.Fr 



In a previous article [Ij we presented two analytical approaches for studying 
the derivative expansion of the exact renormalization group equation (ERGE, 
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for reviews and recent pedagogical introductions see Pl3] ). The two meth- 
ods, based on the commonly used field expansion, were shown to be more 
efficient than the current approaches which implicitly assumed that 

the simple field expansion converges in [0, cx)[ pQ whereas it does not p]. In 
the first method, introduced in the infinite-boundary condition is explicitly 
accounted for via an auxiliary differential equation (ADE) whereas, in the sec- 
ond, the solutions looked for are approximated by generalized hypergeometric 
functions (HFA). Another method very similar to HFA has almost simultane- 
ously been proposed in [8], it looks for the solutions under the form of Fade 
approximants. Both three methods work well (though the ADE method has 
a wider range of application) but they are rather heavy to implement (see 
table [1]). In the present work we show that a simple conformal mapping onto 
the unit disc of a suitably chosen angular sector of the complex plane of the 
independent variable, compactiffes the originally inffnite integration domain 
so as to make the series of the field expansion in the new variable convergent 
on the whole disc of unit radius. 

The paper is organized as follows. In section 1, the principle of the mapping 
method is introduced with the example of the Wegner-Houghton RG fiow 
equation [9] in the local potential approximation (LPA). For the fixed point 
solution of this equation in three dimensions, one approximately knows the 
location of the closest singularity in the complex plane of the independent 
variable We show that the best convergence properties provided by the 
method correspond to the largest angular sector compatible with the analytic- 
ity of the solution in the original variable. The calculations of the eigenvalues 
with the mapping method is shown to be easy and we provide the best esti- 
mates ever obtained up to now of the fixed point solution and the eigenvalues. 
In section 2, we consider the Wilson- Polchinski RG fiow equation [lltil2j in 
the LPA. This equation allows us to illustrate the efficiency of the method for 
different values of the dimension 2 < d < 3. Again, we provide the best results 
ever obtained up to now in three dimensions. We determine the locations of 
the critical and multicritical fixed points for (i = 3, 8/3 and 5/2 together with 
the associated eigenvalues for d = 3 and 8/3 with an excellent accuracy. We 
pursue the determination of the critical fixed point for values of d very close to 
2. (At d = 2, the type of solutions we track disappears.) In section 3, we look 
at the second order of the derivative expansion O (5^) by considering explicitly 
the Morris RG equations [13] in three dimensions. These equations are much 
more difficult to treat than the previous ones, even in the LPA. Nevertheless, 
we are able to determine both the fixed point and the eigenvalues with an 
accuracy approaching that obtained with the shooting method P^fTl] . An es- 
timate of the subcritical "odd-exponent" is obtained for the first time O {d"^) . 
It, however, does not compare favourably with existing estimates [15'J16]. We 
discuss the probable reasons of this disagreement. Finally we summarize and 
conclude. 
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1 The Wegner-Houghton flow equation in the LPA 

As detailed in [Tj, to which article the reader is invited to refer for some ba- 
sical definitions if necessary, the study of the existence of fixed points and of 
their stability in the derivative expansion of an ERGE, amounts to look for 
regular solutions in e [0, cxo[ of coupled nonlinear ordinary differential equa- 
tions (ODE). Here is the (constant) scalar field. Hence the two boundaries 
associated to the ODEs under study are (see [1]): 

(1) (f) = where the symmetry of interest is imposed to the solution 

(2) (p = oo where a specific behaviour in approaching this point is imposed 
to the solution 

This problem of solving differential equations with boundary conditions may 
be numerically studied using a shooting or a relaxation method but these 
methods are not always efficient. Instead, one has often recourse to an an- 
alytical method based on an expansion of the solutions looked for in power 
series of about a certain value (the origin 0q = [4J , or the minimum of the 
potential 05] )• One may also consider 0q as an adjustable parameter [B] with 
a view to improve the convergence of the series. Unfortunately these series do 
not converge in the whole range G [0, cxd[ [6ll0l| and the condition at the 
second boundary cannot be explicitly imposed. 

At this stage, it is useful to consider a concrete example. 



1.1 Fixed point equation 

Let us take as a paradigm the Wegner-Houghton equation in the LPA |9]: 

U = \n[l + U"]+dU + (^l-^^(t)U' , (1) 

in which d is the spatial dimension that we shall set equal to three in this 
section, U stands for the potential U {(f>,t), with t = — In (A/Aq) where A is 
the "running" momentum scale of reference compared to a fixed scale Aq, and 
il = dU/dt, U' = dU/d<i), U" = d'^U/d<p^. 

For the fixed point equation with d = ?>: 

In [1 + U"] + 3f/ - |f/' = , (2) 
the conditions at the boundaries for the solution U (0) are: 
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[/(0)=7, (3) 
f/'(0) = 0, (4) 
f/(0) ^ G0^ (5) 

0— >00 

in which 7 (or G) is the integration constant the value 7* (or G*) of which has 
to be determined. The connection parameter 7* (or G*) is often considered 
as a substitute to the solution U* {(f)) because this latter is deduced from the 
knowledge of 7* (or G*) by a simple numerical integration of the ODE. 

For convenience, instead of 7, we shall deal with r = U" (0) from which, 
according to ([2l [3l Hj), one deduce 7 via the relation: 

7 = -^ln[l + r] . 



Let us expand U {(p) about (/) = up to a finite order 2M. Accounting for the 
conditions (131H1) at the first boundary (the origin = 0), it comes: 



1 1 

Um (0) = -- In [1 + r] + -rz + ^ a„ (r) , (6) 

"J n=2 

z = ^\ (7) 

in which we have introduced z for convenience and the coefficients a„ (r) are 
determined as functions of r so that the fixed point equation ([2]) is satisfied 
order by order in powers of z. For example, the two first coefficients read: 



02 [r 



0-3 (r) 



r(l + r) 

r (1 + r) (1 + 7r) 
360 



It remains to determine the value r* of r so that the condition ([5]) at the 
second boundary is satisfied. To impose this condition, we need an evaluation 
for a large value of z of the Taylor polynomial of degree M introduced in ([6]). 
This evaluation is possible if the radius of convergence Rc of the Taylor series 
in powers of z is infinite. However, in general, in the derivative expansion of 
the ERGE, this is not the case: Rc is always finite due to the presence of 
singularities in the complex plane of z. For the unique nontrivial fixed point 
solution of ([2]), knowing the corresponding r*, the singularity the closest to the 
origin has been numerically located in the complex plane of z at the following 
point Zq [To] : 
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|zo| = 9.7344, (8) 

Arg(;zo) = 0.514 vr = a^TT, (9) 

which imphes that 

/?c = 9.7344. (10) 



From the known locahzation of the singularity the closest to the origin, we 
may assume that the solution U* (0) we are looking for is analytical in the an- 
gular sector which symmetrically straddles the positive real axis, has its vertex 
located on the negative real axis at — l-^ol and its angle equal to Arg (2:0) (see 
figure [1]). This assumption imphes that none of the eventual supplementary 
singularities lies inside this sector. 

The basis of our approach consists in exploiting the analyticity of the solution 
U* looked for in the range z G [0,oo[. Actually, if this solution exists then 
there is an angular sector involving the whole positive real axis of the complex 
plane of z, and including the origin, in which U* is analytic. We thus introduce 
a conformal mapping of an angular sector which symmetrically straddles this 
positive real axis onto the unit disc centered at the origin of the complex plane 
of a new variable w defined by: 



w 



;i + z/R) 



:iii 



[l + z/RY 

in which R and a characterize the position of the vertex and the angle of an 
angular sector of the complex plane of z as shown in figure [TJ 



The inverse transformation is: 



R 



l+w 
1 — w 



(12) 



In particular, if we adjust the angular sector in such a way that the presently 
known closest-to-the-origin singularity [see equations (jH [9])] lies right on one 
of its edge i.e.: 



« = a^ = 0.514, (13) 
R = Rc, (14) 

then the resulting Taylor series of the solution U* expressed in powers of w 
should converge on the whole disc |w| < 1. 

Consequently, to obtain an auxiliary condition for determining the value r*, 
we proceed as follows. 
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Fig. 1. When the interior of an angular sector involving the positive real axis of 
z defines a region of analyticity for the function [/* ((/>) (with z = (f^), it may be 
mapped onto the interior of a circle of unit radius for the complex variable w defined 
by eq (fTTj) . In the present case the singularity zq of U* {(p) the closest to the origin 
determines the analyticity domain provided the other singularities lay outside the 
angular sector containing the origin. 

We deduce from ((61 [T2|) the following Taylor polynomial: 

1 

Um (w) = -- In [1 + r] + ^ Cn (r) , (15) 

n=l 

where Um {w) = Um [(p (w)] with (w) given by ([71 [T2l) . 
Then, from the convergence of the resulting series in powers of w, we may: 
(1) either, impose when w 1 (the counterpart of 2; — > 00), that Um{w) 
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satisfies the counterpart of ([5]), 
(2) or, in a less rigorous but often efficient manner, implicitly exploit the 
convergence of the Taylor series in powers of w on the whole unit disc 
by imposing, as proposed in [TT] , that the M^^ coefficient cm ij) van- 
ishes (this assumes that the remaining higher terms of the series do not 
contribute too much). 

Choosing one of these two variants yields an auxiliary condition under the 
form of a polynomial equation for r the zeros of which are candidates for the 
solution r* we are looking for. 

For the values given by f[TU| [T51 [T^ . and the variant ([T]) of the method as 
example, the real zeros obtained are distributed almost similarly to those 
displayed in figure 1 of [1]. An important difference with [1] however, is that 
it is easier to follow the zero of interest as M grows. The estimation we obtain 
this way for r* is excellent since we get 28 stabilized digits for M = 145: 

r* = -0.4615337201162071199657576484. (16) 



It appears that the values of R and a which correspond to the location of 
the closest known singularity provide the best convergence for r*. Of course, 
this situation is particular because, thanks to the results of [10], we know the 
location of the closest singularity. In general this is not the case and, as we 
shall illustrate in section I2.1.H the only information practically accessible is 
a rough estimation of Rc whereas ac is not known. It is thus justified to try 
different values of R and a. 

Let us define a practical measure of the accuracy of the results as M grows 
by the number of stabilized digits of r*: 



AT,, 



-log 



1 - 



M-l 



M 



(17) 



Figure [2] illustrates well the fact that, if one chooses R^ Rc and a ^ then, 
in general, one gets a worse convergence towards the (hopefully true) value of 

But, there are also other values for the couple {i?, a} which provide essen- 
tially the same convergence, in particular when R is chosen greateiPl than 
Rc- Notice that, in such a case, a must be smaller than Uc otherwise the sin- 
gularity located at zq would have an image in the interior of the unit circle 

^ Taking into account only the criterium of analyticity, one may choose any value of 
R and a such that R < Rc and a < Oc- But, it often appears that the simultaneous 
largest possible values of R and a yield the best convergences. 
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Fig. 2. Study of the fixed point of the Wegner-Houghton equation in the LPA 
{d = 3): approximate number A^'acc of stabihzed digits on [see pT|) ] as func- 
tion of the degree M of the Taylor polynomial in w, for different values of the 
couple {R,a} which characterizes the conformal mapping (jlip . The most efficient 
convergence is clearly obtained for values of R and a corresponding to the location 
of the singularity the closest to the origin [see ([8l[9|)]. 



in the complex plane of w. These favourable cases corresponding to R > Rc 
exist because the closest singularity is not located on the negative real axis 
of z. With the Wilson- Polchinski equation [TTp2] . studied in section ([2]), one 
encounters situations where the condition R < Rcis necessary, indicating that 
the closest singularity is right on (or very close to) the negative real axis of z. 
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1.1.1 Eigenvalues 



A linearization about the solution U* (0) of ([2]), using U {(j),t) = U* (0) + 
e e'^^g (0) with e a small parameter, provides the eigenvalue equation: 

YTU^~2^ +(3 -A) (7 = 0, 

which may be studied by expanding both U* and g in powers of (or z) as 
done above in (P, ED for the fixed point. 

Accounting for the arbitrariness of the normalization of the eigenfunction, the 
conditions at the origin depend on the symmetry of the solution looked for: 

g (0) = 1, g' (0) = , even case , 
g (0) = 0, g' (0) = 1 , odd case , 

whereas the condition at infinity is: 



Once r* is fixed to the value (ITB]) . the coefficients of the series for g depend only 
on A which plays the role of the connection parameter in the preceding fixed 
point equation. We perform the conformal mapping z ^ w defined by (11 II) 
with R and a fixed to the preceding values, namely R = 9.7344 and a = 0.514 
and, because the power of in f ITS]) depends on A we use the variant ([2]) of 
the method, although we could use again the variant ([T]) with the condition 
(6 — 2A) g — (pg' = imposed when w — 1 on the transformed series. The real 
zeros of the resulting polynomial in A are displayed in figure [3] which shows 
an exceptionally clear distribution of the spectrum of the eigenvalues. 

For simplifying the notation, we generically note u the inverse of each posi- 
tive eigenvalue A and u the opposite of each negative one in the even case, 
respectively u and oj in the odd case (eventually with an index when there are 
several eigenvalues of the same kind for each case, see [1] for more details): 

z/ = Y>0 , a; = -A>0. 



Using an expansion up to M = 104, we get an estimate of u with six supple- 
mentary digits compared to the already excellent determination by Aoki et al 

V = 0.68945905616213484062727. 
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Fig. 3. Eigenvalue spectrum (even case) of the Wegner-Houghton equation in the 
LPA (d = 3) hnearized about the Wilson-Fisher fixed point. The open circles are 
the real zeros (as function of the degree of the Taylor polynomial) of the auxiliary 
condition provided by the mapping method (presently using variant ([1])). The hor- 
izontal lines represent the estimates of the eigenvalues (even case). This extremely 
clear picture illustrates well the facility for determining the eigenvalues. 



At order M = 40, the estimates of the eigenvalues for both the even and odd 
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cases are the following: 



V = 0.68945905616, 



0.59524, 



ooi = 1.691340 



2.8385, 



tU2 = 3.999, 



5.19, 



tU3 = 6.40, 



7.7, 



U4 = 9.0. 



The values of the subleading critical exponents are in agreement with the 
previous estimates [18]. As in the case of u, the accuracy of each of them may 
be easily improved by considering higher values of M. 



2 The Wilson-Polchinski flow equation in the LPA for 2 < d < 4 

The Wilson-Polchinski flow equation [111112] in the LPA reads: 



By considering this equation we aim at illustrating other aspects of the method, 
in particular its efficiency compared to other procedures (when d = 3) but also 
the way it works when several fixed points exist (when 2 < c? < 3). A recent 
study of this equation for generic d may be found in [19] and the estimates 
that we provide below compare favourably to those obtained in |20j . 



2. 1 Fixed point equation 
The fixed point ODE reads: 





(19) 



the solutions of which we are looking for have the following propertieJ^: 



^ In [1] the parameter A was noted b. 
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U{0) = k, 
U' (0) = 0, 



+2) ^ o 



,2(d-2)/(d+2) 



(20) 
(21) 

(22) 



where k (or A) has to be fixed to the value of the connection parameter k* 
for A*). 



As in the preceding section, to determine /c*, we consider the following Taylor 
polynomial of degree M in 2; = 0^: 



M 

UM{(p) = k + Y.b„{k) z^, 

n=0 

z = 4>\ 

where the coefficients 6„ (A;) are determined as function of k such that the 
fixed point equation (IT^ is satisfied order by order in powers of z. Then we 
perform the conformal mapping fllip . Since we have a priori no information 
on the analyticity properties of U* as function of z, we consider R and a 
as free parameters. We proceed by trial and error to find their best values: 
we look at the variation of the properties of convergence of the estimate k\,j 
obtained at order M by changing the trial values of R and a. Let us illustrate 
this procedure in the case d = 3 where the value of k* is already known with 
accuracy [21J. 



2.1.1 d = 3 

We use the variant ([1]) adapted to (!22|) by imposing the counterpart of dU* /dz - 
1/2 when w — > 1. Setting arbitrarily R = 1 and a = 1, we easily distinguish, 
among all the real zeros of the polynomial auxiliary condition obtained at 
order M, a generic zero which presumably will converge to k* when M —>■ 00. 
The observed convergence is not excellent yet, but it is sufficient to allow a 
rough estimate of the radius of convergence Rc for the expected k*. Following 
the d'Alembert or Cauchy rules, this estimate may be obtained by looking at: 



Rl,n 



bn 



bn 



+1 



or: 



R2,n — 

for a given value of k and large values of n 



l/n 



In general, we have proceeded that way, by trial and error, to estimate Rf.. In 
the present case, however, the radius of convergence of the fixed point series 
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method 


K 


time 


AUrL; 


u.u r Diyy4uuoizoDo 


iozo.o4 


Fade 


0.07619940081205 


1364.73 


HFA 


0.076199400812340 


138.58 


mapping 


0.0761994008160 


2.00 



Table 1 

Comparison between estimates of the connection parameter k* , of the Wilson- 
Polchinski RG equation in the LPA (d = 3), obtained using different efficient ana- 
lytical methods at order M = 25 of the Taylor polynomial [compare to (j23l) ]. The 
"time" given in the third column is a CPU time (in seconds) corresponding to the 
calculation, on the same computer, using each method. 

is known to be [1]: 

Rc = 5.72167. 

By setting straight out R to this value (to save time), we get a better accuracy 
on the estimation of k* than with R = 1, and on varying a, we observe that 
the accuracy is better and better when a grows. Finally, as shown in figure H] 
the best result is obtained with a = 5/2. This high value of a suggests that 
the singularity is located on the negative real axis and is unique. Indeed when 
we try values of R larger than Rc, the zero of interest disappears contrary 
to the preceding Wegner-Houghton case for which the singularity the closest 
to the origin lies outside the negative real axis. The fact that we can choose 
a > 2 (when R = Rc) without spoiling the property of convergence of the 
mapping method, suggests that this singularity is unique. 

We finally get an accuracy of 49 digits for M = 120: 

k* = 0.07619940081234064145788536913234906280801814336214±6 X 10"^°. 

(23) 

Table [T] shows a comparison of efficiencies at M = 25 between the four meth- 
ods ADE, HFA, Fade and the present mapping method. This latter method 
appears to be, by far, the most efficient. 

The mapping method provides also a global explicit representation of the so- 
lution via the Taylor polynomial in powers of w re-expressed in terms of (the 
HFA and Fade methods give also a global representation of the solution but 
the procedure is much more heavy). From this global representation, knowing 
k*, it is not difficult to estimate the value of the connection parameter A* 
which appears in (1221) . With three more terms in the asymptotic behaviour of 
U* than in (122!) . we have been able to obtain the following result: 

A* ~ -2.3184, 
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Fig. 4. Study of the fixed point of the Wilson-Polchinski equation in the LPA for 
d = 3: approximate number A'^acc of stabihzed digits on k'^j [similar to (fT7|) ] as 
function of the degree M of the Taylor polynomial in w, for different values of a 
and for R = 5.72167. 

which is very close to the value we have determined with the shooting method 

m- 

A* ^ -2.3183. 



Once the fixed point has been determined, the eigenvalues are looked for with 
the same couple {R = Rc, a = 5/2) and, as in the preceding case, we easily get 
the eigenvalues with high accuracy. We may thus give the best estimates ever 
obtained up to now (see [T] for the definitions of both the eigenvalue equa- 
tion and the exponents, and compare the following values with the estimates 
obtained in [21]): 

• for M = 75 in the even case [d = 3): 
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i/ = 0.649561773880648017614299724015827 ±2 x 10" 
u;i = 0.6557459391933387407836879749684 ±2 x 10"^^ 
UJ2 = 3.180006512059167532314140242 , 
^3 = 5.912230612747701026351105 , 
UJ4 = 8.796092825413903643907 , 
UJ5 = 11.798087658336857239 . 

• for M = 69 in the odd case [d = 3): 

ui = 1.8867038380914203710417873172 ± 5.3 x 10"^^ , 
U2 = 4.524390733670772780436353 , 
COS = 7.3376506433543135387526 , 
coi = 10.2839007240259581722 , 
^5 = 13.3361699643459431 . 



2.1.2 d = 8/3 

A new non-trivial fixed point emerges from the Gaussian fixed point when d 
takes on a value below each threshold [22J: 

P 

where p takes on the integer values 1,2,3,-- - ,oo. Each one of these dimensions 
corresponds to a monomial in powers of 0, in the formally expanded potential, 
becoming relevant with respect to the Gaussian fixed point. Hence for 3 < d < 
4, there is only one nontrivial fixed point: the Wilson-Fisher fixed point [23] 
which controls the behaviour of any Ising-like critical point. Below d = 3 
appears a new fixed point which controls the tri-critical behaviour. Down to 
d = 8/3 there are two fixed points. 

To determine these two fixed points, we proceed similarly to the previous case 
d = 3 hj trial and error in order to determine the best values for R and a. 
It appears that we easily determine the critical fixed point and get a rather 
good convergence for R = 3.5 and a = 2. We finally obtain with M = 120 the 
following estimation: 



k*. = 0.16736641293800245119399231271370187±3 X 10 



-35 



It is worth indicating that the radius of convergence of the series for this value 
of k* is: 

i?e = 3.539 . 



An estimation of this connection parameter by a shooting method gives the 
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following result: 



k* = 0.16736641293. 



For the same values of R and a, and the variant ([2]), we get the following esti- 
mates for the (critical) eigenvalues (without trying to optimize the accuracy): 

• in the even case for M = 50: 

u = 0.76520486063609135 ± 1.8 x 10^^^ , 
ui = 0.81489673946448431495 ± 7.4 x 10"^^ , 
UJ2 = 3.32048207983 , 
COS = 6.045184748 , 
cui = 8.9242446 , 
a;5 = 11.92109, 

• in the odd case for M = 60: 

ui = 2.032987010721739017 ± 5 x 10"^^ , 
tU2 = 4.66036178907299, 
0^3 = 7.46810111086, 
C04 = 10.409630969 , 
C05 = 13.4559703 . 

The tri-critical fixed point is more difficult to determine than the critical one 
because a clear identification of the correct zero is only possible above M ~ 20, 
a value which is greater than that observed in the case of the critical fixed 
point. Nevertheless, we get a good estimate for R = 19.5 et a = 0.47 (variant 
([2D, and M = 120): 

k; = -0.0152088617493 ± 9.2 x 10"^^ , 
whereas by the shooting method we get: 

k; = -0.015208861395 , 



The radius of convergence of the series in this case is: 

R^ ~ 19.4. 



For the same values of R and a, and the variant [2l we get the following 
estimates for the (tri-critical) eigenvalues (without trying to optimize the ac- 
curacy) : 

• in the even case: 
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z/1 = 0.501489558 , 

U2 = 1.0479506 , 
6^1 = 0.343226, 
6^2 = 1.83594, 
6^3 = 3.484, 
uj4^ = 5.26 , 
UJ5 = 7.2 , 

• in the odd case: 

i>i = 0.659916, 
6^1 = 1.0682, 
6^2 = 2.643, 
6^3 = 4.36, 
UJ4 = 6.27 , 

6^5 = 8.4. 

As expected, compared to the critical fixed point, the tri-critical fixed point 
has one additional positive even eigenvalue (A = 1/ z/2), and also one additional 
positive odd eigenvalue (A = l/i^i). 

2.1.3 Lower dimensions 

The calculations are more and more difficult as d decreases. This is because the 
radius of convergence Rc (d) for a given kind of fixed point (critical or multi- 
critical) also decreases. Consequently, in approaching d = 2, the determination 
of the critical fixed point itself is made more and more difficult. In fact, at 
d = 2, the regular solution we are searching (which satisfies the condition 
( [22I) when (j) — > 00) disappears. Indeed, the same reasoning as that followed 
by Morris in [21], shows that for d = 2 the only solutions to the fixed point 
equation in the LPA are either singular at finite or periodic (see also [25]). 

For the sake of shortness, we have simply tested that we can determine the 
critical fixed point (we have not considered the eigenvalues despite the fact 
that they are easy to determine) down to dimensions very close to d = 2. 
Table [2] shows some of the results for k*. It shows also that the radius of 
convergence of the original series in powers of decreases when the location 
of the minimum of the critical fixed point potential increases. Consequently, 
because this property is presumably general in the derivative expansion, the 
method of [5,6J based on an expansion about this minimum is presumably 
doomed to failure before reaching d = 2. 

Regarding the fixed points of higher criticality, we have not tried to determine 
them below d = 5/2. For this value of d we get the locations of the fixed points 
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1.41 


21/10 


0.96514877904597 


0.75 


2 


2.214026 


1 


202/100 


1.8411524065675 


0.55 


2 


2.370782 


0.77 



Table 2 

Determination of the connection parameter k* of the critical-fixed-point solutions of 
([T9|) for decreasing values of the spatial dimension d. <pQ stands for the location of the 
minimum of the fixed point potential (C/*'(0o) — 0) ^ determined from the present 
mapping method. rc{(j)) = \l Rdz) is the radius of convergence of the original series 
in powers of (j). Notice that rc{4') decreases when the location of the minimum of 
U* increases and becomes smaller than \(I)q\- R and a are the parameters of the 
conformal mapping defined in (jlip . 

displayed in table [31 

We show in figure [5] the global solutions we get from the Taylor polynomials 
in w for the fixed points when d = 5/2. They have the right n-well potential 
form with n growing with the number of directions of instability of the fixed 
point in agreement with [22] . 



3 The Morris equations (LPA and second order of the derivative 
expansion) 

At second order of the derivative expansion [O (9^) in short], in addition to the 
potential U (0, t), a new function appears: the coefficient Z (0, t) of the kinetic 
term of the action S [0]. The flow equations for U and Z are coupled and the 
reparameterization invariance of the ERGE [26] is broken [27I1T3] together 
with the so-called scheme invariance (invariance with respect to a change of 
the cut-off function which introduces the running momentum scale of reference 
A). This is, at least, the case for the Wilson- Polchinski ERGE. Actually the 
two invariances are not independent [2S] and if, instead of the action S [0] , 
one considers the Legendre transformed or effective action F [(f], then the 
breaking of the invariances may be reduced to a single one [2] . One may even go 
further and find an appropriate cut-off function so that the reparameterization 
invariance is effectively restored. Morris [13] has thus obtained the following 
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type of fixed point 


1 mapping 
1 shooting 


min(M) 




(Rc 


a 


critical 


1 


0.2529954929769066647431 
0.252995492980 


15 


< 


I 2.71 
[ 2.71 


2 


tri-critical 


< 


-0.043027 
-0.043026994781 


50 


< 


9 

~ 13 


0.6 


quadri-critical 


1 


0.0030001 
0.003000121150 


40 


1 


30 
~ 27 


0.5 



Table 3 

Determination of the values of the connection parameter A;* for the three non- 



trivial fixed points of (I19p when d = 5/2. min(M) indicates the minimal value 
of the degree M of the Taylor polynomial at which, in the mapping method, we 
can unambiguously identify a zero that one may easily track when M grows. The 
estimates obtained by using the shooting method are shown on the second line 
within a bracket. R and a are the parameters of the conformal mapping (lllh . the 
Rc are estimated values of the radius of convergence of the solution in terms of the 
original independent variable z = (f)'^. 

two equations (written below for d = 3 and keeping the original notation J^: 



V 



ir=(l-V4) 



_^^JI^= + 3V -Ul + ri)xV' 
kJv" + 2v^ ^ 



1 2AKK" - 19 {K') 



(24) 



48 



if 3/2 (^^,, ^ 



3/2 



1 bW'KWK + 57 {K'Y + {y"'Y K 
5 {y"'f K + 2V"'K'VK + {K'f 



12 



1 



K iy" + 2v^' 



7/2 



- (1 + ri)xK' -rjK , 
2 



(25) 



Going from S [cp] to T [if] changes the couple {U, Z} into {V, K} 
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Fig. 5. The mapping method provides a precise global explicit information on the so- 
lutions looked for. Here are shown the n-well potentials as obtained from the Taylor 
polynomials in w {(j)) of the three fixed-point solutions found for the Wilson-Polchin- 
ski RG equation in the LPA when d = 5/2. 

where x stands for and r/ is a parameter to be adjusted in order to find a 
fixed point. 

Let us first consider the determination of the Wilson-Fisher fixed point. 
The fixed point equations correspond to: 



V^ = 0, (26) 
i^ = 0, (27) 

where V and K stand for the rhs of fl24l [25 
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We look for solutions under the form of even functions of x which must satisfy 
the following conditions: 



V{0) = k, 


(28) 


V (0) = 0, 


(29) 


KiO) = Ko, 


(30) 


K'iO) = 0, 


(31) 


V (x) ~ Gix^ , 


(32) 


K (x) ~ G2X i+" . 


(33) 



and which have no singularity in the whole range x G [0, 00 [. 

The two parameters k and rj (or Gi and G2) have to be adjusted to the two 
connection parameters k* and 77* respectively (or G^ and G2) which define 
the fixed point solution of intereslQ- By virtue of the reparameterization in- 
variance, Kq may be set arbitrarily equal to one without changing the value 
T]* (and of the eigenvalues associated to the linearization of the flow equation 
about the fixed point). 

As in the case of the LPA, we express the two functions as Taylor polynomials 
of degree M in z = x'^: 

M 

l^M(x)=fc+^6« (A;,r/)^", 

n=0 
M 

n=0 

z = x'^ , 

in which the coefficients b^^ {k, t]) and b^^ {k, r]) are determined so as to satisfy 
the fixed point equations fl26l [271) . To determine the values k* and rj* using 
these two polynomials, we need two auxiliary conditions that we deduce from 
the conditions fl32| [331) and the analyticity of the two functions V (x) and 

K [x) in an angular sector of the complex plane of the independent variable 

2 

z = X . 

We thus introduce the same conformal mapping z ^ w as defined by ([TTj) and 
we adopt one of the two variants defined in section [TTTl to get the two auxiliary 
conditions either by directly imposing on the Taylor polynomials the counter- 
part of the asymptotic conditions (1521 [551) when w 1 or, less rigorously, 

^ The parameters k* and rj* are related to the parameter a and tj used in [13] by 
rf = T] and k* = ' 
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by assuming the simultaneous vanishing of the two coefficients ll^l} {k, rj) and 



3.1 LP A 



Let us ffist study the LPA of fl2^ 1251) which then reduce to the following 
unique equation: 

V = — -i + W--V' . 

^JV" + 2 2 

The numerical study of this flow equation is much more complicated than that 
of the Wegner- Houghton or Wilson-Polchinski cases. The largest degree of the 
Taylor polynomial reasonnably accessible is, in general, smaller than in these 
preceding cases. Nevertheless, using the mapping method, we have correctly 
estimated the fixed point location and the eigenvalues with a better accuracy 
than the results available in the litterature [TSfl^ (obtained by the shooting 
method, see table H]). For a better comparison, we have improved those latter 
estimates (see table Hj). 

For the determination of the fixed point in the LPA [d = 3), the maximum 
degree of the Taylor polynomial that we have employed is M = 82. The best 
estimation of k* , given in table HJ is obtained with R = 2 and a = 1/2. The 
radius of convergence Rc of the original series in powers of z is roughly equal 
to 2. These facts suggest that the singularity the closest to the origin lies close 
to the imaginary axis in the complex plane of z. 

The process of determining the eigenvalues is similar to the cases already 
described in the preceding sections. 

In the LPA, the eigenvalues have been estimated with shorter series than the 
fixed point, M = 40 and M = 60 for the even and odd cases, respectively, 
whereas R and a remained unchanged. The clear distribution of the eigenval- 
ues shown in figure [3] is preserved at the LPA and O (<9^). This particularity of 
the mapping method is very interesting since it is not always easy to clearly 
determine the spectrum of the eigenvalues (using the shooting method for ex- 
ample, see table Hj). In the LPA the mapping method provides essentially the 
same accuracy than the shooting method. 



3.2 Second order 



At the next order 0((9^), the main difference with LPA is the number of 
equations: two instead of one. For the sake of shortness we shall not describe 
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k* 



UJl 



0J2 



UJl 



UJ2 



LPA 
mapping method 
shooting (this work) 
shooting (previous) 



0.2753644064810282 
0.2753644064810124 
0.275364406 



0.660389431 

0.660389431331 

0.660389 

0.6285575 

0.6285575035 

0.6285 

3.04801 
3.048005033 



1.8124863608 
1.812486361 



4.32251050 
4.3225104975 



0(92) 
mapping method 
shooting (this work) 
shooting (previous) 



0.258216 
0.258214391 
0.258 2144 

0.053941 

0.0539320839 

0.05393208 

0.618063 

0.61806 

0.6181 

0.8964 

0.897227 

0.8975 

1.71 
1.7 



0.86562 
0.865569 



2.8 



Table 4 

Estimates of the critical and subcritical exponents of the Morris equations (j24l I25p 
at d = 3 in the LPA, and at the seconc?§rder of the derivative expansion. On the 
right of each bracket are given the results obtained by the mapping (first line) and 
the shooting (second line) methods used in this work and, when they exist, previous 



again the different steps of the calculations. 

For the fixed point, the maximum degree of the Taylor polynomial we have 
attained is M = 17. Moreover, for small values of M, it has been difficult to 
locate the correct solution of the fixed point equation because the auxiliary 
polynomial conditions have a large number of real zeros. Nevertheless we have 
been able to distinguish the correct zeros. The values of k* and rj* displayed 
in table H] have been obtained with R = 2.5 and a = 2. We have estimated 
the common radius of convergence to be Rc {z) ^ 2.5. For other values of 
R and a we have observed different kinds of convergence towards a unique 
couple of values. So, we are able to provide a rough estimate of the error on 
the evaluation of the O (9^)-fixed point by the mapping method. For M = 17, 
it comes: 

r = 0.258204 ±0.000023, 
r/* = 0.05388 ±0.00011. 

Regarding the eigenvalues, the largest degree of the Taylor polynomials con- 
sidered is M = 25. The parameters R and a were not necessarily fixed to the 
values mentioned above for the fixed point. The eigenvalues displayed in table 
Hlhave been obtained with R = 2 and a = 2. In addition to those eigenvalues, 
we have observed in the even case the presence of a very small eigenvalue 
1.5 X 10~^). It corresponds to the expected zero eigenvalue associated with 
the reparameterization invariance. Its absolute value may be seen as a rough 
measure of the accuracy of the calculations. 

It is worth indicating that our O (9^)-results are better than the estimates: 



obtained in |6j using an expansion about an adjustable non-zero value of the 
field. 



3.3 Remark 

Let us emphasize that table H] displays for the first time, for the Morris equa- 
tions fl24f25l) . values of the eigenvalues in the odd case {uJi, 0J2) together with 
the first estimate of UJ2. The exponent ui is of particular interest since it is 
actually the first time that it is evaluated at the second order of the deriva- 
tive expansion of an ERGE. Moreover, this quantity is not well determined 



77 = 0.05425, 
z/ = 0.617476, 
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since only two estimations are available: from the Wilson ERGE [TT] using the 
scaling field method in [T5] : 

uji = 2.4 ±0.4, 

and from the e-expansion up to the order in \l6l from which we have ex- 
tracted the value: 

= 2.34 ±0.49. 

These two results are close to each other and clearly compatible. This is not 
the case of our present evaluation from the Morris equations (B^ Indeed, 
assuming that the derivative expansion converges, we may roughly give an 
error bar on the estimate of ui. Using the two values of table HJ it comes: 

cji = 1.34 ±0.5, 

a result which is incompatible with the two preceding estimates. It is thus 
likely that the present O (9^)-calculations do not provide the correct value of 
uji, either because the derivative expansion itself does not converge or sim- 
ply because the Morris equations are not adequate to determine correctly the 
values of the critical exponents (at least at low orders of the derivative ex- 
pansion). We think that the second possibility is the right explanation. This 
is because one has already observed in the results of [29] a bad convergence 
on uJi and also on rj* as functions of the number N of components of the field 
in the O (iV)-symmetry whereas another choice of cut-off function, as done 
in [30], seems to be better adapted (see [31j for example). It is very likely 
that the fundamental reason why the Morris equation is not efficient must be 
looked for in the constraint imposed to the O (c}^)-equations with a view to 
satisfy explicitly the reparameterization invariance. Perhaps this constraint is 
too strong at this order. 



4 Summary and conclusion 

We have presented a new efficient analytical method for solving two point 
boundary value problems when one of the two boundaries is located at infinity. 
The method is based on the analyticity of the solution looked for in an angular 
sector which contains the positive real axis including the origin. This angular 
sector is then conformally mapped onto the unit disc so that an expansion in 
powers of the new independent variable w of the solution looked for yields a 
convergent series on the whole disc. The conditions at the infinite boundary 
may then be safely imposed on an approximate solution written as a Taylor 
polynomial in w estimated when w ^ 1. 
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We have illustrated and tested the method on three different kinds of ODEs 
that one encounters in the derivative expansion of the ERGE: the Wegner- 
Houghton and the Wilson-Polchinski equations in the LPA and, at second 
order of the derivative expansion, the Morris equations. For the first and sec- 
ond examples (LPA) we provide the best estimates of both the fixed point 
connection parameter and of the critical exponents ever produced up to now 
from the ODEs considered. For the second example we are able to follow the 
fixed point solution from the dimension d = 3 down to 2 where the kind of 
solutions we track disappears. For the intermediate dimensions d — 8/3 and 
d — 5/2, we are able to locate the existing multi-critical points and the critical 
exponents attached to them with great accuracy. The difficulties increase with 
the multiplicity of the fixed point and with the decreasing of d. The study of 
the Morris equations is more demanding since it involves two ODEs (second 
order of the derivative expansion) which are particularly difficult. We show 
that the method works again well. We are able to estimate subcritical expo- 
nents that were not calculated previously using the derivative expansion. The 
estimation of one of these subcritical exponents does not compare favourably 
with values coming from the e-expansion and the scaling field method for 
studying the ERGEs. We give some possible explanations why the two equa- 
tions considered are presumably not the best choice for calculating the critical 
exponents. 

Gomparisons with other analytical methods show that the mapping method 
is actually very efficient. It may be useful in the study of more complicated 
systems of equations obtained from the derivative expansion. 

Obviously, it can also be employed efficiently to study analogous kinds of 
ODEs that appear in other fields when one of the boundaries is at infinity. 
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